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The weak localization correction and the mesoscopic fluc- 
tuations of the polarizability and the capacitance of a small 
disordered sample are studied systematically in 2D and 3D 
geometries. While the grand canonical ensemble calculation 
gives the positive magnetopolarizability, in the canonical en- 
semble (appropriate for isolated samples) the sign of the effect 
is reversed. The magnitude of mesoscopic fluctuations for a 
single sample exceeds considerably the value of the weak lo- 
calization correction. 

PACS numbers:73.20.Fz,73.23.-b,73.61.-r 



I. INTRODUCTION 

The phenomena of weak localization (WL) and meso- 
scopic fluctuations in disordered systems have been inten- 
sively studied during the last fifteen years, mainly in con- 
nection with transport properties of these systems . 
For these phenomena the role of the electron-electron 
interaction is just in setting the length scale (phase 
breaking length), below which the electron wave func- 
tion preserves its phase coherence. Systems with a size 
L less than are called mesoscopic systems. 

In this paper, we consider quantum corrections and 
mesoscopic fluctuations of the two other characteristics of 
a mesoscopic system, where the electron-electron interac- 
tion is essential: polarizability and capacitance. The for- 
mer quantity can be measured by putting a sample into 
a capacitor, while the latter one determines the charging 
energy which shows up in the Coulomb blockade experi- 
ments Q. 

The flrst quantum calculation of the polarizability of 
a small metallic particle was obtained in a seminal paper 
by Gor'kov and Ehashberg (GE) [Q. It rehed on the fol- 
lowing two assumptions concerning statistical properties 
of energy levels and eigenfunctions in disordered systems: 

1. The single-particle energy spectrum exhibits the 
same statistics as the eigenvalue spectrum of ran- 
dom matrices from the Gaussian Ensemble of the 
corresponding symmetry; 

2. Exact single-particle eigenfunctions tpkif) and 
ipiir), which are close enough in energy, are cor- 
related in the following way: 



{rk{r)Mr)Mr'W{r'))^^^ = UD{r,r'). (1) 
Here the average is defined as 

mr)Mr)Mr'Mir')),,^^ (2) 
{J2rkir)Mr)Mr'W{r')S{e - e,)<5(e + u; - q)) 

k^l 

k^l 

V is the sample volume, and the diffusion propagator Hd 
is a solution to the diffusion equation. 



DV'^TlD{r,r') = (Triy)"^ [5{r - r') ~ V~ 



(3) 



with the boundary conditions Vnllr) = 0. The first of 
these conjectures was proven by Efetov and the sec- 
ond by the authors More precisely, it was shown 
in 1^ that for the energy difference much less than the 
Thouless energy, lu E^, 



{ru{r)Ur)Mr'Wi{r')), 
= kd{r - r') +ni)(r,r'). 



(4) 



where kd{r) = (7rzy)~^ (ImG^(r))^ is a short-range func- 
tion {G^{r) being the retarded Green's function) explic- 
itly given by 



fcd(r) = exp(-r//) 



Jgippr), 2D 



(5) 



The short-range part kd{r — r') of the correlation func- 
tion (^ was not taken into account by Gor'kov and 
Eliashbcrg, but it would give only a small correction to 
their result . 

Based on these conjectures, GE concluded that the po- 
larizability for very low frequencies a; ^ A (A = (i'V)~^ 
is the mean single-particle level spacing; v is the density 
of states per spin) is enhanced in comparison with the 
classical value ao ~ the enhancement factor being of 
order (kR)^, where 



(87re2i.)i/2^ 3D 
ine^iy, 2D 



(6) 



is the inverse screening radius. Although the original 
paper [Q gave a new insight into the field (called later 
mesoscopic physics), and had a substantial impact on 
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the further development of the condensed matter physics, 
this result for the polarizability is incorrect for the follow- 
ing reason. The paper by Gor'kov and Eliashberg does 
not take into account the effects of screening: they calcu- 
late the polarizability in response to the local field rather 
than to the external one As was found in Refs. [p|Jlo| 
(see also the screening restores the classical value 

of the polarizability, thus reducing the quantum effects 
to a relatively small correction. Evaluation of this cor- 
rection was recently attempted by Efetov , who com- 
bined the non-perturbative calculation of the polarization 
function with the electron-electron interactions taken 
into account in the RPA approximation. Since the value 
of the quantum correction depends on the presence (or 
absence) of the time reversal symmetry, it was denoted 
by Efetov as a "weak localization correction to polariz- 
ability" ; we are following this terminology in the present 
article. However, he estimated incorrectly the contribu- 
tion of the short-range term in Eq.(^, which made him 
to conclude that the weak localization correction is dom- 
inated by this term. As we show below, this is not the 
case if the system size exceeds considerably the mean 
free path. More recently, Noat, Reulet, and Bouchiat 
p3[ presented a perturbative calculation of the weak lo- 
calization correction to the polarizability in a particular 
geometry of a narrow 2D ring. They considered both 
canonical (CE) and grand-canonical (GCE) ensembles, 
and concluded that the correction to the polarizability 
is parametrically suppressed in the CE. While essentially 
confirming their GCE result, we disagree with the above 
statement concerning the CE. We show below that the 
effect in the CE is of the same (up to a coefficient of order 
one) magnitude than the GCE one, but has an opposite 
sign. 

As was realized by Berkovits and Altshuler fluc- 
tuations in the polarization function lead to mesoscopic 
fluctuations of the polarizability of the sample. They 
considered a specific thin film geometry and identified 
the four-diffuson diagrams giving the leading contribu- 
tion to the fluctuations. We will follow their approach 
when studying the polarizability fluctuations in 2D and 
3D geometries. 

Along with the polarizability, we consider another 
quantity characterizing a mesoscopic system, the capaci- 
tance C. It determines the charging energy e^/C, which 
manifests itself in the I — V characteristics of a quan- 
tum dot in the Coulomb blockade regime. In particular, 
the charging energy represents the main contribution to 
the threshold voltage in the excitation spectra and to 
the distance between adjacent conductance peaks in the 
addition spectra. Statistical properties of the Coulomb 
blockade I — V characteristics are attracting a great re- 
search interest now JTs)"^, which motivated us to con- 
sider the WL correction and the mesoscopic fluctuations 
of the charging energy. In addition, the capacitance de- 
termines the low-frequency behavior of the impedance of 
mesoscopic systems [ p2| , p3[ . 

Therefore, the purpose of the present paper is to study 



systematically the WL effects and mesoscopic fluctua- 
tions of the polarizability and capacitance in 2D and 3D 
geometry. Where it is necessary, we refine results of pre- 
vious research. We show that the polarizability and the 
capacitance can be treated on the same physical grounds. 
We will also find a simple relation between the magni- 
tude of the WL correction and that of the mesoscopic 
fluctuations. The electron-electron interaction is taken 
into account in the RPA approximation, which works for 
K < pp, pp being the Fermi momentum. We consider the 
case of low temperature T <C A (thus setting T = in all 
formulas), and study both grand-canonical and canonical 
ensembles. 



II. WEAK LOCALIZATION CORRECTION TO 
THE POLARIZABILITY OF SMALL PARTICLES 

We consider an isolated disordered metallic particle 
(3D or 2D) placed into a uniform external frequency- 
dependent electric field E{uj). We assume that the sys- 
tem is diffusive, I <^ L, where / and L are a mean free 
path and a typical size of a particle, respectively. In the 
RPA approximation the potential distribution $(r) and 
the electron density p{r) in the particle obey the Poisson 
equation (e — ~\e\ being the electron charge), 

= { i othe'iwis^ ' 
in combination with the equation 

p{r)^-2e I n{r,r')^r')dr'. (8) 
Jn 

In 2D we use the following convention throughout the 
paper: r = {x,y) denotes the coordinates in plain, and 
z the transverse coordinate. The Laplacian A is always 
a three-dimensional operator, A = A3 = A2 + d"^. Fur- 
thermore, n is the polarization function (per spin), which 
can be readily expressed through the Matsubara Green's 
functions, 

n(r , r',u;) = -Tj2 (G(r , r', ie„, + iuj.^) 

(9) 

or, in terms of the retarded and advanced Green's func- 
tions G^'^ {r , r' , e) , 

Uir,r',Lj) 

= S/ de{np{e){G^{r,r',e + u;)G''{r',r,e))- 

- np{e + uj){G^{r, r' , e + cj)G^(r', r, e)) + 

+ [np{e + u;)- nj.(e)](G^(r, r', e + u;)G^{r', r, e))} 

= n^^(r, r', + n^-^(r, r', + n^^(r, r', lu). (10) 
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While and H"^^ can be calculated in the usual impu- 
rity perturbation theory, the evaluation of 11^"* for low 
frequencies, w ^ A, requires a non-perturbative treat- 
ment. This was done in Ref. in the framework of the 
supersymmctric sigma-modcl approach, and we present 
only the results here: 



n 







- I delm{G"{r,r',e)f 

7^ J-oo 



+ — [(RcG^(r,r',0))2 - (ImG«(r, r', O))^]; 

ZTTl 

U^\ry,u) ^-y,-^ {(ReG«(r,r',0))^ 
V Zm 

+S{uj) (ImG«(r,r',0)}2 + {7riy)^{l + 5H)nz5(r, r')} • 

(11) 



Here S'(w) is a correlation function of the zero- 
dimensional sigma-model, S{llj) — —{QllQb^) (in nota- 
tions of Rcf. PI ) , given explicitly by 

(12) 



S{lj) = 1 H — exp 



TTILU 



sm ■ 



It is related to the two-level correlation function R2{oj) 
as R2{ijj) — Re[l + 5(w)]/2. Now we decompose the 
polarization function into frequency-independent and 
frequency-dependent parts: 



n(r, r', u:) = no(r, r') + ni(r, r', oj), 
no(r,r')--Im / de(G^(r, r', e)}^ - 

ni(r,r',c.) = -;^(l + 5H) 

■vfXlDir^r')] 



V 



(13) 



2m 

[(ImG^(r,r',0)}2 
^A{u^)[Ur - r') 



ni5(r,r') 



(14) 



where we defined A{Ld) = [i-KLu /2/S.){l + S{Ld)) and intro- 
duced the function kd{r — r') = (ImG^(r, r' , 0))^ given 



exphcitly by Eq.(|l). Note that the formulas du]), ([Til) 
( p^ ) are written for the case of the unitary ensemble 
(broken time-reversal symmetry due to the presence of 
a strong enough magnetic field), which was considered in 
Ref. . Generalization to the orthogonal ensemble (un- 
broken time-reversal symmetry) is straightforward and 
results in the following modification of the factor A{lo) 
in Eq.dH) for Hi: 



^(^) = ^(l + 5(c.))-l. 



(15) 



where now 

S{u) = 1 ■ 



2ie** sins 



2i — 

ds 



f°° e*"** 



and s ~ TTw/A. In the limit of low frequency, a; <C A, 
the factor A{uj) is equal to 



A(.«A,.-i.{;J; , ,17, 

where [3 is the usual parameter equal to 1 (2) for the 
orthogonal (resp. unitary) ensemble. 

Now we turn to the calculation of the dipole moment. 
The general expressions were obtained in Rcf. jl^]; we 
present here the derivation for the sake of complete- 
ness. We consider the frequency-dependent part Hi as 
a perturbation and expand functions p(r), ^(r), and the 
dipole moment, 



d ~ e J rp{r)dr, 



with respect to Hi. In the zero-order approximation we 
obtain 



do = -2e rno(r,r')$o(r-')rfrdr', 
Jn 



(18) 



where the potential <i>o(^') satisfies Eqs. (0), (^) with 11 
replaced by IIq. 



%(r) = 




FIG. 1. The electrostatic potential and the electron den- 
sity in the RPA approximation. The external dashed line is 
the bare potential of the electric field — Er, the loops with 
and 1 denote the contributions Do and IIi to the polar- 
ization function, respectively. The wavy line is the Coulomb 
interaction. 

It is easy to check that the first term in the expression 
( [Tsl ) for no(r, r') gives ly after integration over one of 
the coordinates, the integral being dominated by the dis- 
tances of order of the Fermi wave length, \r — r'\ ^ Xp- 
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Assuming the screening length (which sets the scale for 
the variation of the potential ^{r)) to be larger than Xp, 
we can thus replace this term by a delta-function: 



no(r,r') ^i^[S{r-r')-V-^] 



(19) 



This approximation for the polarization function Hq leads 
to the simple relation between the potential ^Q{r) and 
the excess density po{r) (Thomas- Fermi approximation) 



%{r) -(2et.)-Vo(r-), 



(20) 



and consequently to the following equation for the po- 
tential <I'o('"), 



A$o(r-) = 



_ r K,^<^>o{r)9nir) , 



3D 
2D 



(21) 



with K defined in Eq.(^ and the boundary condition 
$0(7- ^ 00) = Er + const. In Eqs.(|2^), we have 
chosen the arbitrary additive constant in the definition 
of $0 in such a way that ^o{'r')dr = 0. 

In the following we consider particular geometries of a 
3D spherical sample of a radius R (to be referred as 3D), a 
circle of a radius R in the in-plane electric field (2D), and 
a quasi-two-dimensional sample of a thickness h {0 < z < 
h) and an area S in the field directed along the z-axis, i.e. 
transverse to the sample (Q2D) The polarizability tensor 
aij is generally defined as d{uj) — a(uj)E(Lu). Eqs. ( |l^ ) 
and (|l|) yield the classical polarizability ||2^ equal in the 
limit kR <C 1, fv/i ^ 1 to (for all the geometries under 
consideration the dipole moment is directed along the 
field, and the tensor is reduced to a scalar). 



3D 

ao ^ <, Sh/An, Q2D 
(4/37r)i?3, 2D 



(22) 



Now we turn to the corrections due to the function Hi . 
One obtains (Fig.l) 



diE = 2e^ / drdr'^a{r)Ili{r,r')^oir'), 



(23) 



in fuU accordance with Ref. [|2|. Note that Eq. (|2|) 
depends explicitly on the symmetry of the system with 
respect to the time reversal, and therefore constitutes the 
WL correction to the polarizability. As follows from Eq. 
(p^), this correction consists of two contributions. The 
first one (to be referred as ais) is due to the short-range 
contribution to the polarization function (the first term 
in the brackets, Eq. (|lj)), and the second one (ai/j) is 
due to the diffusion contribution Hd- 

To evaluate the second term, we use the expansion of 
Hd in the eigenfunctions of the Laplace operator (j^aif) 
with the boundary conditions V^q — and the corre- 
sponding eigenvalues — Eq, 



n£,(r,r') = {-kDv) 



:-Va(r-)0„(r'). (24) 



For the purposes of estimate, one can use the following 



expressions valid for I <^ \r 
of the system), 



' r'l ^ L (L is a typical size 



{2TT^vD)-^\TL[L/\r -r'\], 2D 
{ATT'^vD\r ~r'\)-^, 3D 



(25) 



Now we evaluate and compare both contributions 015 
and aiD in 3D and 2D systems. Since the structure of 
the potential <I>o is different in 3D and 2D cases, these 
should be treated separately. 



A. 3D geometry 

For any 3D geometry with nL ^ 1 the expression for 
<l'o('") can be written in the form 



*o(t") = — </3(r||)exp(-Kr_L), 

K 



(26) 



with ip being some function of magnitude unity. We have 
introduced a transverse coordinate r± (r for the sphere, 
z in the case of a disk), and the vector ry of coordinates 
along the surface of the sample. Note that according to 



Vif\\) = -47recro(r||)/£; , 



(27) 



where e(To(T'||) is the charge density on a surface of an 
ideal conductor induced by the electric field E. It can be 
found by the methods of the classical electrostatics . 
In the integral for a is, 



"IS 



2e'^iyA{uj) 



j d^rd^r'^^{r)ka{\r-r'\)Mr'), 



(28) 



both points r and r' lie in fact in the layer of thickness 
<Si I near the surface of the sample. One can then 
integrate over the transverse coordinates and reduce the 
remaining double surface integral to the integral over one 
coordinate only. We obtain 



1 



1 



A{u;)\n{Kl) / d'n^\n) 



^A(.)ln(.0. 



(29) 



In Ref. 1 12 1 the kernel in Eq. ( pq ) was incorrectly re- 
placed by a 5-function, which led to an overestimate of 
the contribution ais by a factor of nl{\n{nl))~^ » 1. 
On the other hand, for the term due to the diffusion. 



aiD 



2e^v 
VE'^' 
x<i>o(r'). 



A{uj) / drdr'$o(r)nz5(r,r') 



(30) 
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we obtain, using the estimate (E 



1 



aiD 



vDk- 



:A{uj) ~ ais 



L 



l\n{niy 



(31) 



We see that in a diffusive system of a size L ^ l\n{Kl), 
the diffusion contribution auj dominates, in contrast to 
the conciusion of Ref. At the same time, if the sam- 
ple size L is comparable to the mean free path I (which 
happens e.g. in ballistic systems with surface scatter- 
ing), the short-range contribution is parametrically of the 
same order (in fact, even larger by a logarithmic factor) 
as the diffusive one. As expected the WL cor- 

rection is small in comparison with the classical polariz- 
ability ao, 



and the diffusion term dominates for h > l\n{Kl). In 
the opposite case of a thin sample the short-range con- 
tribution is the leading one. In particular, for h < I we 
find 



ai 



ais = A{uj) 



(38) 



As is seen from the above formulas, the relative magni- 
tude of the weak localization correction is rather low for 
both 3D and quasi-2D (with the field direction normal to 
the plane) geometries, so that the experimental observa- 
tion of the effect in these cases may be problematic. The 
effect is much more pronounced in the 2D case, which we 
consider below. 



ai/ao 



g{KL)' 



(32) 



g ~ 2'iTvDL being the dimensionless conductance. The 
calculation of the numerical coefficient for the WL cor- 
rection to the polarizability requires the exact expansion 
(^) . For the particular spherical geometry the potential 
$0 has a form 

3-£/-R 

$o(r) = exp(-K(i?-r))cos6', R-r<^R, 

TTKr 

(33) 



and we obtain 



ai = 



1.36 



(34) 



According to Eq.(|17|), the WL correction to the polar- 
izability is negative. The value of the correction in the 
presence of strong magnetic field (unitary symmetry) is 
smaller (twice as small for zero frequency) as without 
the latter (orthogonal symmetry). The experimentally 
measured magnetopolarizability as , defined as 

as =a(S) -a(0), (35) 

is therefore positive, in agreement with Ref. [T^ . 

B. Quasi-2D geometry (transverse field). 

We consider now a quasi-two-dimensional sample of a 
thickness h ^ and an area S ^ with the elec- 
tric field directed transverse to the sample plane. Then 
Eq.(p6|) for the potential reduces to 



$o(r) = -(-e-^^-t-e- 

K 



i{h-z) 



(36) 



If the sample is relatively thick, h ^ I, the same consid- 
eration as for the case of a spherical shape yields 



ais = A{uj) 



Iii{kI) 



C. 2D geometry (in-plane field). 

In contrast to the 3D case, the potential $o in the case 
of a 2D sample in the in-plane electric field is a smooth 
function of coordinates, with the characteristic scale set 
by the sample size R. Therefore the kernel in the integral 
( p8| ) , which has a support of order I , can be replaced by 
a ^- function, 

{pFr)~^ exp{—r/l) « 2'Klpp^S{r). 

This gives an estimate 

ais ~ l{pFKy^A{uj)\n{R/l). 

On the other hand, for the diffusive term (|3^) the esti- 
mate (^5|) implies 

aiD ~ RHKg)-^A{cu) ^ aisiR/lfiHR/lJr' ■ 

Similarly to the 3D case, the diffusion term ais domi- 
nates for R ^ I. The relative magnitude of the quantum 
correction can thus be estimated as 

ai/ao l/gKR, 

with g — 2t:i'D — kpl/2. 

For the particular case of a circular geometry, the po- 
tential (f>o is given in the polar coordinates (r, 6) by |23| 



$o(r) = -2E{TTK)-'^rcos9{R^ - r^) 



-1/2 



(39) 



An exact calculation gives the value of the quantum cor- 
rection 

ai{uj) = 1.53R^{kpfI)-^A{uj), 
and the relative magnitude of the correction is 

1 



ai/ao — 3.6 x 



kPfIR 



Aico) 



In a recent paper, Noat, Reulet and Bouchiat (NRB) 
|l3t proposed a geometry of a narrow 2D ring (radius 
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R, width W ^ R) as more favorable for observation of 
the effect. In the in-plane electric field the ring becomes 
polarized with the one-dimensional (i.e. integrated over 
the ring cross-section) charge density 



ER 



■ cos 9, 



e\n{R/W) 
and the classical polarizability given by 



ln(i?/t^) ■ 



Calculating the quantum correction, we find again that 
for a diffusive ring i? 3> the contribution aio domi- 
nates and gives 



i?4 



ai 



vDW^k\t?{R/W) 
The relative magnitude of the correction is 



Oil/oiQ 



1 



TTgWK\u{RlW) 



A(^), 



(40) 



(41) 



where g is now the quasi-one-dimensional conductance 
g — vDW/ R. These results for the ring geometry are by 
and large in agreement with those found by NRB [ p^ . 
Actually, NRB express the polarizability in terms of the 
exact eigenfunctions of electrons, conceptually similarly 
to the original GE calculation, and then perform the im- 
purity averaging using the semi-classical expression for 
the correlation of the exact single-particle eigenfunctions 
(see e.g. p],p^). This calculation yields correct results 
for the following reasons. First, the short-ranged terms 
in Eq. (^, omitted in this calculation, turn out to be 
unimportant for the WL correction to the polarizabil- 
ity. Then, the exact expression for the long-ranged (dif- 
fusive) part of the eigenfunction correlator (|l|) coincides 
with the semi-classical result even for w <C A, where the 
latter generally is not expected to be true. This has been 
proved and discussed previously by the authors in Ref. 

®-. . 

Similarly, we can consider a quasi-one-dimensional 
strip of width W and length L ^ W oriented along the 
electric field direction (which we choose to be the z-axis) . 
Again, the sample polarization is described by the one- 
dimensional charge density 



Ez 



eln{L/W) ' 
yielding the classical polarizability 

" 12\niL/W) ■ 
The quantum correction is now equal to 



ai ~ aiD 



SOKWgln^iL/Lu) 



(42) 



where g = 2'kvDW/ L is the dimensionless conductance. 
We obtain 



27r 



1 



5 KWg\^{L/W) 



(43) 



Thus, we have found that the WL correction to the 
polarizability can be quite appreciable in 2D (circle) and 
especially in quasi-lD (ring or strip) geometries, which 
gives a possibility of its experimental observation. These 
conclusions are in full agreement with those of NRB, Ref. 



D. Canonical ensemble 

The results obtained above were derived for the grand 
canonical ensemble, where the chemical potential is fixed 
by an external reservoir. In the Appendix A we present 
the calculations for the canonical ensemble, which is more 
appropriate for the problem in question [p6| . We show 
(Eq. (p2[)) that the CE magnetopolarizability differs from 
the GCE result by the coefficient —2.75. This means that 
although the magnitude of the CE effect is the same as 
in the GCE, the sign is opposite in the CE case: the 
magnetic field suppresses the polarizability. 



III. WEAK LOCALIZATION EFFECTS IN THE 
CAPACITANCE 

A. Definitions 

A natural definition of the capacitance in an open sys- 
tem is 



= edQ/dfi, 



(44) 



where Q is the total charge of the system. This capac- 
itance determining the low-frequency transport proper- 
ties of the system was studied thoroughly by Biittiker 
and coworkers in Refs. |2^,0 (where it was called "elec- 
trochemical capacitance"). An explicit calculation p^ ] 
yields 



e 

a, 



1 



Vv{Ef)' 



(45) 



where Cg is the geometrical capacitance, determined 
from the equations of the classical electrostatics |2^ with 
corrections due to the screening effects, and iy{Ep) is the 
density of states at the Fermi level. The average of the 
second term in the right-hand side of Eq. (|4^) is the 
mean level spacing A, which is usually much less than 
the average of the first term. However, the fluctuations 
of the second term are important p3t . 

Being formally applied to the closed system, Eq. (|4^ ) 
yields infinite fiuctuations of the charging energy. Indeed, 
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in Ref. [Q, where an attempt to calculate the fluctua- 
tions of the compressibility of a closed system has been 
made, the integral over energies diverged, and the au- 
thors could get a flnite result only by cutting it off at 
energies of order A. 

A proper generalization of the definition (Q) for a 
closed system is its "discrete" version. 



(46) 



where fJ-{N) is the chemical potential of a closed system of 
N electrons. The quantity (^6|) has an important physical 
meaning: it is equal to the spacing between two consec- 
utive peaks in the addition spectrum of a quantum dot 
in the Coulomb blockade regime. Statistical properties 
of these spacings were studied experimentally in Refs. 
[|6|-|l8| and theoretically in Ref. 

Similarly to the case of an open system (Eq.(|45|)), the 
peak spacing (46) can be decomposed into two parts: 
level spacing A^r and (usually much larger) contribu- 
tion associated with Coulomb interaction (denoted Ei 
in Ref. [^). The main contribution to the latter (and 
thus to the Coulomb blockade peak spacing) is given by 
the charging energy Ec defined as a constant part 
of the effective two-particle interaction potential U (r, r') 
of the electrons in the sample: 



Ec = eye = V 



I drdr'U{r,r') 
Jn 



(47) 



In particular, the WL correction to the charging energy 
calculated below is the dominating term for the WL in 
the Coulomb blockade peak spacing, and can be, in prin- 
ciple, measured as the magnetic field dependence of the 
peak spacing. In contrast, fluctuations of the charging 
energy do not give the dominant contribution to the fluc- 
tuations of the peak spacing, see Sections IV, V, and Ref. 




FIG. 2. RPA approximation for the two-particle potential 
U. The potential Uq is given by the same sequence of diagrams 
as the potential $o in Fig. 1 provided the external dashed line 
is replaced by the Coulomb interaction (the wavy line). 



B. Weak localization effects 

The potential U{r, r') can be found in the RPA ap- 
proximation (see Fig. 2). As in the case of the polar- 
izability, it is convenient to split the polarization function 
into two parts (^9|) and (|l^). Assuming the low frequency 
limit a> <t; A, we replace the function A{u;) by its zero- 
frequency value, Eq. (p7|) . In the zcroth order in Hi one 
should solve an equation 

Uo{r,r') - Voir -r')~2 f dr^dr2Vo{r - n) 

X no(ri,r2)[/o(r2,r'), 1^0(0 = A- (48) 

This equation was solved for an arbitrary closed system 
in Ref. H]. The result is 

C/o(r, r') = U + Mr) + $o(r-') + C/«(r, r'). (49) 

Here U = (e^/C)o is a constant, corresponding to the 
charging energy calculated in the Thomas-Fermi approx- 
imation; Uk is the usual screened Coulomb potential, 
shifted by a constant so that J drUt^{r,r') = 0, while 
$0 is the contribution due to the excess positive charge, 
moved towards the boundary of the system after an extra 
electron is added to the system. For the sphere (3D) and 
circle (2D) geometries this potential has an explicit form 



~ , , f e2(Ki?2)-iexp(-K(i?-r)), ?,D 

$o(r) = const-|^,j2^^5_,^^,^_^,y,/,^^^ 2i^ ' 

where the constant is chosen in such a way that 
/ drl>o(r) = 0. 

In the first order in Hi we obtain (Fig. 2) 

{/i(r,r') = 2 y"dridr2{/o(r,ri)ni(ri,r2){/o(r2,r'), 

and, taking into account that the integral of Hi over any 
of the coordinates is zero, we write the corresponding 
contribution to the charging energy (which constitutes 
the WL correction) in the form 



2 / dridr2l'o(r-i)ni(ri,r2)$o(r-2). (50) 



Due to the structure of the function Hi (Eq. (||)) there 
are two contributions to the WL correction: one comes 
from the short-ranged term, and another one is related 
to the diffusion. A comparison of these contributions 
can be carried out exactly in the same way as was done 
for the polarizability, and it turns out that the diffusive 
term dominates if the sample size exceeds considerably 
the mean free path. Thus, we obtain finally 



-^j = I dridr2<^>oiri)TlD{ri,r2)<i>o{r2). 



(51) 
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The calculation for the particular geometries of a sphere 
(3D) and disk (2D) gives 



-1.32/3-V-l(pi.i^)-^ 
-0.010/3-V-l(pF^?)"^ 



3D 
2D 



(52) 



The small coefficient in front of the 2D expression is an 
artifact of the specific circle geometry. 

The weak localization correction suppresses the charg- 
ing energy, i.e. enhances the capacitance. The magnetic 
field suppresses the capacitance. Both in 2D and 3D 
cases, the WL correction to the charging energy can be 
estimated as 

e2\ A 



(53) 



The WL correction can be in principle extracted from 
the measurements of the magnetic field dependence of 
the capacitance, though its rather small value may make 
such a measurement problematic. 

As for the polarizability, the CE calculation yields an 
additional factor -2.75 in Eqs. d^) and (H). Thus, 
in CE the magnetic field suppresses the charging energy, 
enhancing the capacitance. 



IV. MESOSCOPIC FLUCTUATIONS 

Mesoscopic fluctuations of the polarizability and the 
capacitance can be calculated in a similar way. In addi- 
tion to the average polarization function (p^), ( p^ ) there 
is also a random part 11^ (r,r') with the zero average, 
giving rise to the fluctuations of these quantities. Since 
the integral of 11^ over each coordinate is zero, we imme- 
diately arrive to the expressions for the random parts of 
the static [ p8| polarizability Ur and the charging energy 
{e^ /C)r in precisely the same way as Eqs. ( p3| ) and (Im 
were obtained: 



r 

ar = -^ / drdr'<^o{r)Iir{r,r')<^o{r') 



and 



2 / drdr'^Q{r)I{r{r,r')^Q{r'). 



(54) 



(55) 





a) 



b) 



FIG. 3. Diagrams for the fluctuations of the polarization 
function. The double dashed lines denote the diffusion prop- 
agators. The counting factors are 2 (a) and 4 (b) for the 
unitary symmetry, and 4 (a) and 8 (b) for the orthogonal 



Thus, the mesoscopic fluctuations of both quantities 
are determined by the fluctuations of the polarizability. 
To calculate them, we perform the perturbative calcu- 
lation. As we will see, in the case of fluctuations the 
whole range of energies A < e < contributes, and this 
fact justifies the perturbative calculation (in contrast to 
the weak localization correction, which is determine d by 
low energies). Following Berkovits and Altshuler p4[ , 
we identify the four-diffusion diagrams (Fig. 3) as giving 
the main contribution to fluctuations of the polarization 
function. 

We obtain 

(n,(ri,r3)n,(r2,r4)) = -^^Re / tdeD,{r^,r2) 

xD^{r2,r3)D^{rs,r^)D^{ri,r:i). (56) 

Here 12//3 is a combinatorial factor (a number of 4- 
diffuson diagrams); the function is given by the fol- 
lowing expression. 



1 

1^ 



TivTVoiry), A<e<E, , (57) 



and decreases for e > Ec- If the integral over e in Eq. 
( p6| ) diverges, the cut-off at energies of order of the mean 
level spacing A, where the perturbative expression ( [s^ ) 
ceases to be valid, should be introduced. 

A naive estimate suggests that the leading contribution 
to the fluctuations of the polarizability as well as of the 
capacitance can be found by substituting the zero-mode 
contribution (the first term in rhs of Eq. (^)) for all 
the four functions D^. However, since the integrals of 
the potentials $o and $0 are zero, this term vanishes. 
For typical geometries the main contribution is given by 
the term where two of the functions are replaced by 
the zero-mode result, while two others are represented by 
the diffusion propagator (the second term in rhs of Eq. 
(p7|)). Cutting off the logarithmically divergent integral 
at A from below and at the Thouless energy Ec from 
above, we obtain 



/ 2\ 12e' 21 



J dridr2$o(n)nz3(ri,r2)$o(r2) 



(58) 



and 



12 2, 



V 



dr 1 dra $0 (n) Br, (ri,r2) $0(^-2) 



(59) 
2 



Comparing these results with the expression for the WL 
corrections, Eqs. ( ^0| ) and (|5l|), we obtain a general for- 
mula relating WL corrections (calculated within GCE) 
to the mesoscopic fluctuations of the same quantity. 
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(a) = (3/? In .g) 



and 



(3/3 In g) 



ai\ 



(60) 



(61) 



Note that in contrast to the WL corrections, the meso- 
scopic fluctuations are not expected to be sensitive to the 
GCE/CE difference, since the integral ( |5^ ) is determined 
equally by all energies A < e ^ £^c- 
For particular geometries we obtain 



r.m.s a 



and 



C 



/3 



1/2 



1/2 



4.7iri(pi.K)-2, 

5.30i?2(Kp^Z)-\ 



3D 
2D 



3D 
2D 



The remaining terms in Eq. (|56|), with all the four 
functions replaced by diffusion propagator, can also 
be easily estimated. Their contribution is of the same or- 
der as Eqs. (|60|), (|6l|), but without the logarithmic factor 
in the numerator. Thus, for the r.m.s. of the charging 
energy in addition to the term of order ^ A{\ngy^'^ /g 
(two zero modes) we obtain a correction of order ~ A/g 
(no zero modes). 

Amazingly, the relations between the weak localiza- 
tion correction and the amplitude of mesoscopic fluctua- 
tions, Eqs. (|60| ) and (|6ll), have a universal form, the same 
for 2D and 3D systems. We should note, however, that 
these results are not applicable for the case of polarizabil- 
ity fluctuations of a quasi-two-dimensional sample of the 
thicknes s h ^ I and area S ^ in the transverse field 
(see Se cJiI Bh calculated previously by Berkovits and Alt- 
shuler [|14||. In this case, the contribution from the terms 
with two zero-modes, Eq. ( pSh, w hich can be easily cal- 
culated with the use of Eq. ([36|), yields the r.m.s value 
of the polarizability fluctuations. 



r.m.s(Q;^2)) — 



(3/?- 1 Ing) 1/2 h 



ng 



where g = 2'KvDh. On the other hand, the contribu- 
tion where all the four functions are replaced by the 
diffusion propagators is Q 



8/3-i7r3/2 51/2 



i.e it is larger by a factor ^ S'i/2/[/i(lng)i/2], and for this 
particular geometry represents the leading contribution 
to the fluctuations of the polarizability. 



V. CONCLUSIONS 

In this paper, we have calculated the weak localization 
(WL) correction to the polarizability and the capacitance 



of a disordered sample and the mesoscopic fluctuations of 
these quantities. The WL corrections originate from the 
GrGa term in the polarization function, which depends 
on the presence or absence of the time-reversal symme- 
try. A change of the polarization function influences the 
screening and, consequently, the polarizability and the 
capacitance. We find that in the grand canonical ensem- 
ble, switching on the magnetic field leads to a positive 
correction to the polarizability and negative one to the 
capacitance. In the canonical ensemble the magnitude 
of the effect is the same (up to a numerical coefficient 
^ 2.75), however the sign is reversed. 

Calculating the mesoscopic fluctuations of the polar- 
izability, we find that for typical geometries they are re- 
lated to the value ai of the WL correction as follows (see 
Eq. d^)): 

r.m.s.(a) = (3/?lng)i/2|ai|. 

The same conclusion is valid for the capacitance, see 
Eq.(|l]). Therefore, the magnitude of fiuctuations ex- 
ceeds considerably the value of the WL correction. This 
should be contrasted with the relation of the correspond- 
ing quantities for the case of the conductance: 

r.m.s.(g) 1; 

1 , quasi- 1-D 

|gi| ^ < ln(L/0 , 2D 



L/l 



3D 



so that |gi| 3> r.m.s. (g) in 2D and 3D and |gi| ^ r.m.s. (g) 
in the quasi- ID geometry. As our results show, an ex- 
perimental observation of the magnetopolarizability of 
mesoscopic samples requires an experimental setup with 
large number of such samples, which would reduce the 
fluctuations. 

Mesoscopic fluctuations of the charging energy con- 
tribute to the fiuctuations of the conductance peak spac- 
ings in the addition spectra of quantum dots in the 
Coulomb blockade regime [p"5|-p^. However, as follows 
from Eqs.(|5|), (|6l|), the magnitude of these fiuctuations 
is much smaller (by a factor {In g)^/^ /g) than the level 
spacing A. Therefore, the contribution of the charging 
energy fluctuations to the fiuctuations of the peak spac- 
ings is parametrically smaller than the effect of electron 
level fluctuations, which is given by the random matrix 
theory and is of order of A. The charging energy fiuctua- 
tions represent one (but not the only one) of the contribu- 
tions to the enhancement of the peak spacing fiuctuations 
as compared to the random matrix theory. This problem 
was considered in detail in Ref. |20[ . 

Finally, we would like to mention once more that we 
assumed the screening length to be much larger than the 
wave length, or in other words, that Ts ^ 1, where Ts = 
e^/evp- In the opposite case, > 1 (but still below 
the Wigner crystallization threshold), one can roughly 
estimate the result assuming the screening length to be 
approximately given by the distance between electrons. 
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This leads to an enhancement of the above results for the 
WL correction and the r.m.s. amplitude of fluctuations 
by a factor of order of Ts (resp. r^) for the polarizability 
and capacitance, respectively. 
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APPENDIX A. MAGNETOPOLARIZABILITY IN 
THE CANONICAL ENSEMBLE 

In this Appendix we calculate the magnetopolarizabil- 
ity as defined by Eq. ( ^5| ) in the canonical ensemble. 
For this purpose we first rewrite our derivation of the 
Sec. II in terms of exact eigenfunctions and energy eigen- 
values, in the same manner it was done originally by GE 
Q, and later by NRB ||l^. Using the results, derived 
for the correlation of eigenfunctions in Refs. (uni- 
tary ensemble) and Appendix B (orthogonal ensemble), 
we first obtain the weak localization correction in the 
grand canonical ensemble. It is in full agreement with 
the results obtained in Sec. II. Then we generalize the 
derivation to the CE case. 

Our conclusion is that the CE magnetopolarizability, 
ag^, can be easily obtained from the GCE value, ^ 
as follows: 



A ^GCE 



A 



CE 



2.75. 



(62) 



We should stress here that the derivation given below 
uses only the linear response formalism and the prop- 
erties of the eigenfunction and eigenvalue statistics in 
disordered systems. It relies on the fact that that the 
WL correction to the polarizability is determined by the 
energy range where the level correlation is important. It 
can be repeated for the case of the capacitance, where one 
obtains a relation analogous to Eq. (|6^. In contrast to 
this, the fluctuations of both quantities are determined 
by the energies A ^ e <C i?c, where the correlation of 
levels does not play any role, and the difference between 
CE and GCE values is not expected. 



Magnetopolarizability in GCE: derivation a la 
Gor'kov and Eliashberg 

We start from Eq. (p3|), which can be rewritten in 
terms of the exact single-particle states as follows, 

asiuj) ^ ^ dridr2^o{ri)SIl{ri,r2)^o{r2). (63) 



Here we introduced the symbol S denoting the difference 
between the quantities in unitary and orthogonal ensem- 
bles. 

Si- ..) = (.. ■)gue — (■ • ■)goe, 
and the polarization function 11 is expressed as follows, 

m {,'^2 ) 



- (ni + e„ + iO ■ 



(64) 



m and n being the exact single-particle states. Thus, for 
w A we obtain an expression valid both in GCE and 
CE, 

"b(-) = ^^/ E ^^l('J'o)_rV (65) 



In GCE the position of the Fermi level can be arbitrary, 
and we replace the sum in Eq. (65) by an integral with 
the level correlation function i?2(e): 



erfei?2(e)(...) 



(66) 



e„<CF<em 

Using the sum rule for the eigenfunctions. 



we obtain for the magnetopolarizability 
2e2 



as 



(67) 



The above derivation of Eq. (67) is essentially equiva- 
lent to that of Refs. ||,|3|. Now, using Eq. (|79) for the 
case of the orthogonal symmetry and Ref. [29| for the 
unitary one, we write (r = |ri — r2|) 

^'(|V'm(ri)Vm(r2)P). 

" [1 + 2nD{ri,r2)] [1 + 2kdir)] , GOE 
[l + ni3(ri,r2)][l + fcd(r)], GUE ' 

Separating the dominating diffusion term, we obtain for 
the magnetopolarizability in the grand canonical ensem- 
ble 



2e^ 



as 



V^E^A 



J dridr2$o(ri)nz3(ri,r2)$o(r2), (68) 



which coincides with Eq. (30) 
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Canonical ensemble 

To realize the canonical ensemble, we apply the 
method previously developed in Refs. js^js^. Namely, 
we fix the number of electrons to be integer in each in- 
dividual sample, but allow it to fluctuate slightly from 
sample to sample. This type of ensemble is realized by 
pinning the Fermi-level to one of the single-particle levels 

^k- CF = Efc + 0. 

Now instead of Eg. (pq ), one should split the sum over 
energy levels in Eq. (pq) into two. The first contribution 
consists of the terms with m — k, and can be transformed 
to the integral with the use of the two-level correlation 
function i?2 . The rest of the sum requires the three-level 
correlator _R3(0, e, ei), which corresponds to the proba- 
bility to find three levels with energies 0, — ei, and t — ei 
(counted from the Fermi-surface). Thus, we obtain 



/>OC 

E (•■■)=A-W &i?2(e)(...) 

1*00 /•€ 

-f / de / deii?3(0,ei,e)(. 
Jo Jo 



(69) 



and for the magnetopolarizability 



as 



de 



(70) 



i?2(e)A+ / deii?3(0,e,ei; 
Jo 



Since the integral over e converges for e ~ w, we can re- 
place the matrix element by its low- frequency limit. For 
e <C -Ec we obtain for the correlation of the wavefunctions 
with m ^ n (see Eq. (^2|) for the orthogonal symmetry 
and Ref. for the unitary one): 

' kd{r) + [l + kd{r)]nD{r^,r^), GOE 
kd{r)+JlD{ri,r2), GUE ' 

In the cases when the diffusion dominates one can 
replace this eigenfunction correlator for either symmetry 
by n^i. As a result, taking into account Eq. (|68|), we 
obtain Eq. (^), with the coefficient Ace expressed in 
terms of the level correlation functions, 



A 



CE 



-1-2 



ds 



(71) 



i?2(s)+ / dsi(i?3(0,Si,s)-i?2(s)) 




where we have made a change of variables s — 7re/A. 
Note that Eq. ( [7l|) differs from the similar expression 
derived by NRB [ [13[ by the second term in the rhs. 

In the leading approximation the correlation function 
i?2 and i?3 may be taken from the random matrix theory 



P3[ . The first integral in the rhs of Eq. 
calculated analytically and is given by 



1) can be 



7o s Jo s 



-0.367, 



in contrast to the statement of NRB that it is equal to 
— 1/2. We have defined the functions 



sm s 



9{s) = -f{s)- his) 



/(si)dsi 



The second term in the rhs of Eq. ( |71| ) after a lengthy 
algebra can be expressed as follows: 



h= — dsiS{R3{Q,si,s)~R, 
Jo s Jo 



is)) 



ds 

— I dsi{g{si)h{si) - f{s)g{si)h{s - si) 



+gis)fis,)h{s - si) + h{s)fisi)gis - s,)} . 

Calculating this numerically, we find I2 = —1.509 and 
thus. Ace = -1 - 2(/i + h) = 2.753. 

Note that in the above derivation we neglected the con- 
tribution of the so-called Debye processes (relaxation to 
the instantaneous equilibrium distribution due to cou- 
pling with phonons or other possible inelastic processes) 
P4| , p5 13 1 . These processes do not exist in a closed sam- 
ple in the limit of zero temperature (T A) that we are 
considering p6|. 



APPENDIX B. CORRELATIONS OF 
EIGENFUNCTIONS IN DISORDERED 
SYSTEMS: ORTHOGONAL ENSEMBLE 

In this Appendix we derive the expressions for the cor- 
relations of the eigenfunctions in the orthogonal ensemble 
in the same way as these were obtained in Ref. [|9) for 
the unitary ensemble. We restrict ourselves to the terms 
of order g"^. 

Following j2^] we define the eigenfunctions correlators 
(see Eq. d)), 

il{i'i,r2,e) = {\'4'k{ri_)ipk{r2)?) ^ 
/3(ri, rs, e, uj) = (|?/'fc(ri)V';(r-2)P), ^ , k^l, 



(72) 



and 

7(ri,r2,e,w) = (V'fc(n)'0i(n)'0fc('^2)V'r('^2))e,^ , k^l. 
The quantities 77 and (3 are related as follows, 

B = v{ri,r2, e)A-^5{uj) + /3(ri, ra, e, lj)A-^R2{lj) 
= v^ + (27r2)-iRe {(G«(ri, n, e)G^(r2, r2, e + uj)) 
- (G^(ri,ri,e)}(G^(r2,r2,e + a;))}; (73) 
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here the two-level correlation function, 

R2{^) = A2 J(e - ek)S{e + lu ~ ei)\ , (74) 
is introduced. 

The right-hand side of the expression ( [73| ) can be di- 
rectly calculated with the use of the supersymmetry tech- 
nique. For the case of preserved time-reversal symmetry 
(orthogonal ensemble) one obtains 

B(ri,r2,e,c^) = -(27r2)-iRe {(g^i ,i(ri, n) 
X 5fci%i('^2,r-2) + .g6i%i(ri,r2).ggi\bi(r2,ri))^ 

- (5bMi(ri,ri))^ (36i%i(r2,r2))^} . (75) 

Here {■ ■ ■)f denotes the averaging with the action of the 
supermatrix sigma-model 



(•••)f = j DQ[...)eM~F[Q]), 
F[Q] = j dr SiY[D{VQf + 2i{ 



w + iO)AQ], (76) 



where D is the diffusion coefficient, Q = T^^AT is 
a 8x8 supermatrix, A = diag(l, 1, 1, 1, — 1, — 1, — 1, — 1), 
and T belongs to the supercoset space U{2, 2|4)/L/(2|2) x 
U{2\2). The symbol Str denotes the supertrace (trace 
over bosonic degrees of freedom minus that over fermionic 
ones). The upper matrix indices correspond to the 
retarded-advanced decomposition, while the lower indices 
denote the boson-fermion one (here we need only the in- 
dex 61, which denotes one of two bosonic components of 
a supervector). The Green's function g in Eq. (|75| ) is the 
solution to the matrix equation: 



-i(e +--Ho)- -{uj + iO)A + Q/2t 
= (5(r-r'). 



g{r,r') 



(77) 



Expressing these functions through the matrices Q and 
taking into account Eq. (^3|), we arrive at the following 
equation valid in for an arbitrary diffusive system: 



27:' 



Siu;) 



/3(ri,r2,e,w) 



i?2(w) 



:(^^.)2Re{l-(Qii,i(ri)g22,i(r2))F 
-2fcrf(r)(Qii%i(ri)Q2i^,i(ri))p}, 



(78) 



here the function kd is defined in Eq. (||), and r = |ri — 
r2 1 . The separation of the the rhs of Eq. ( [78|) into the 
singular (proportional to S{uj)) and regular parts allows 
one to obtain the quantities a{ri,r2) and f}{ri,r2,uj). 

For the case of a metallic system in the weak lo- 
calization regime, the sigma-model correlation functions 

{QllM(^i)Qllbiir2))F and {QHMQll,Ar2))F can 
be calculated for relatively low frequencies uj ^ Ec with 



the use of a general method developed in Refs. [g7 38 
which allows one to take into account spatial variations 
of the field Q. The results are obtained in the form of 
an expansion in g' 
obtain 



Up to the terms of order g ^ 



we 



{QiiMiri)QitMir2)) 

= 1 - 2R{u) 



4iA 



and 



tt{lj + iO) 



{QllM{ri)Ql\A^2))F = -2 



iA 



-2 Riuj) + 



nD(T'l,T2) 

zA 

Tr(uj + iO 

nz)(ri,r2). 



Tr[uj + iO)^ 

Here the diffusion propagator Hd is defined by Eq. (|^), 
and we have introduced the function R{u!) = [H-S'(w)]/2, 
where S{uj) is given by Eq. (|l|). Note that the two-level 
correlation function, i?2(a^), is the real part of R{uj). 

Now, separating regular and singular parts in rhs of 
Eq. (|7^), we obtain the following result for the autocor- 
relations of the same cigcnfunction, 



= [l + 2fcd(r)] [l + 2Hz5(ri,r2)] 



(79) 



and for the correlation of amplitudes of two different 
eigenfunctions {k ^ I) 

V^{\Mri)Mr2)\'')e,u - 1 = 2kd{r)nDiri,r2). (80) 



The result (79) for ri = r2 is the inverse participation 
ratio previously obtained in Rcf. while that for an 
arbitrary spatial separation was found in the zero-mode 
approximation [g — oo) in Ref. ]39[ |. 

Now we turn to the correlation function 7. Similarly 
to Ref. [E9|, we obtain a relation 



2tt 



Siu;) 



7(^i,'f2,e,w) 



R2{lo) 



A ' A2 

= ~{7Ti.)'Re{{Qll,,{rMlbii^2))F + kd{r) (81) 

X [{Qllir^)QlUri))F + {Qllbiiri)QllMir2))F - l] } • 

Separating again the rhs into the regular and singular 
parts, we recover Eq. (|79| ) and obtain 

V^rk{ri)Mri)Mr2)^!ir2))e.^ 

= kd{r) + [l + kd{r)]nD{ri,r2), k^l. (82) 
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